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ABSTRACT 

Using RXTE/ASM archival data, we investigate the behaviour of the spectral hardness ratio 
as a function of X-ray luminosity in a sample of six transient X-ray pulsars (EXO 2030+375, 
GX 304-1, 4U 0115+63, V 0332+63, A 0535+26 and MXB 0656-072). In all sources we find 
that the spectral hardness ratio defined as F 5 _i 2 kev/-F'i. 33 - 3 kev increases with the ASM flux 
(1.33-12 keV) at low luminosities and then saturates or even slightly decreases above some 
critical X-ray luminosity falling into the range ~ (3 - 7) x 10 37 erg s 1 . Two-dimensional 
structure of accretion columns in the radiation-diffusion limit is calculated for two possi¬ 
ble geometries (filled and hollow cylinder) for mass accretion rates M ranging from 10 17 to 
1.2xl0 18 g s 1 . The observed spectral behaviour in the transient X-ray pulsars with increasing 
M can be reproduced by a Compton saturated sidewall emission from optically thick magne¬ 
tized accretion columns with taking into account the emission reflected from the neutron star 
atmosphere. At M above some critical value M cr ~ (6 - 8) x 10 17 g s 1 , the hight of the col¬ 
umn becomes such that the contribution of the reflected component to the total emission starts 
decreasing, which leads to the saturation and even slight decrease of the spectral hardness. 
Hollow-cylinder columns have a smaller height than the filled-cylinder ones, and the contri¬ 
bution of the reflected component in the total emission does not virtually change with M (and 
hence the hardness of the continuum monotonically increases) up to higher mass accretion 
rates than M cr for the filled columns. 
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1 INTRODUCTION 

Accretion of matter onto a gravitating centre is a ubiquitous phys¬ 
ical process in astrophysics. The matter gravitationally captured 
from the surrounding rarefied medium by a body of mass M (or, 
if in a binary system, escaping from the secondary star as a stellar 
wind or due to Roche lobe overflow) eventually reaches the sur¬ 
face of the body with the velocity of the order of the parabolic one, 
v = y/lGM/R, where R is the radius of the body. This velocity 
becomes very high in the case of accretion onto a compact star. 
For neutron stars with radius R NS ~ 10 km the velocity reaches 
~ 10 10 cm s _1 , so that the accretion energy release is L ~ O.lMc 2 , 
where M is the mass accretion rate and c is the velocity of light. 
These basic facts required the calculation of the accreting mat¬ 
ter braking near the compact star surface. In the pioneer work by 
[ZeTdovich & Shakura] |1969j >, the structure of the shock wave pro¬ 
duced by spherically symmetric accreting matter near the neutron 
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star surface without magnetic field was calculated. Braking of the 
accreting matter was assumed to be primarily due to Coulomb inter¬ 
actions of the infalling protons with the surface, which is justified 
at low accretion rates when interaction of the charged particles with 
radiation is insignificant. 

Immediately after the discovery of X-ray pulsars in 1971 @3 
iacconi et a l.|1971j > it was realized that in bright pulsars the radiation 
plays a crucial role in braking of the accreting matter onto the sur¬ 
face of a neutron star with strong magnetic field. Indeed, in this 
case the matter is canalized by the neutron star magnetic field onto 
the magnetic polar caps, which are characterized by a small radius 
>'o ~ Rns / Ra ~ 10 5 cm (for a purely dipole field), where R A 

is the magnetosphere radius. 

It is now well recognized that once the radiation density above 
the polar cap starts playing the role in the accreting matter dy¬ 
namics, an optically thick accretion column above the polar cap 
is formed (Basko & Sunya ev] 1 976} (see |Mushtukov et al. 1 1 }20 1 5] > 
for the latest investigation). The characteristic height of the column 
increases with accretion rate, most of the emission escapes through 
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Figure 1 . The ratio of the fluxes in 5—12keV and 1.33-3keV ranges (referred to as “hardness ratio” in text) measured with RXTE/ASM for six transient 
accreting pulsars as a function of the total ASM flux in the 1.33-12 keV range. 


sidewalls, and it can be expected that above some X-ray luminosity 
the dependence of the observed properties of the continuum emis¬ 
sion on the X-ray luminosity can become different from that in the 
low-luminosity regime. 


Here we stress that the effects of the strong magnetic field are 
important for the observed properties of X-ray continuum at en¬ 
ergies below the cyclotron line energy due to scattering processes 
which have different cross-sections for different photon polariza¬ 
tion modes. At small accretion rates, the accretion mound can be 
viewed as a thin slab, and mostly ordinary photons (for which 
the electron scattering cross-section is strongly angular dependent 
but frequency independent) escape along the magnetic field lines 
forming a pencil-beam X-ray emission diagram (Basko & Sunyaev 
|1975| >. In contrast, when the accretion column is formed at higher 
X-ray luminosities, mostly extraordinary photons (for which the 
electron scattering cross-section is angle independent but strongly 
frequency dependent) escape sidewall to form a fan-like X-ray 
emission beam ( |Lyubarskii|1986) . 


Presently, there is a wealth of observational data on X-ray pul¬ 
sars that can be used to check different regimes of accretion onto 
magnetized neutron stars. Different types of spectrum-luminosity 
dependence were first discovered based on the behavior of cy¬ 
clotron resonant scattering features (CRSFs) (e.g.|Tsygankov et al. 


2006); 


2012 ) ; 


Stau bert et al.| ( |2007) ; |Klochkov et aT7|4201 l| >; |Becker et al. 


Nishimura 


( (2014) ). 


However, the mechanism of the CRSF generation in accre¬ 
tion columns of accreting neutron stars is still not well understood 
due to a very complicated physics of radiation transfer in super¬ 
strong magnetic fields of neutron stars. For example, recently a 
novel model for CRSF production in spectra of bright X-ray pul¬ 
sars due to reflection of the radiation generated in the accretion 
column from the neutron star surface was proposed (Poutanen et al. 


|2013j l. This model can successively explain the CRSF energy de¬ 
pendence on luminosity observed in outbursts of the transient X-ray 
pulsar V 0332+53 (see also ( |Lutovinov et al.|2015) for further ap¬ 
plications of the model). From the observational point of view, the 
study of CRSFs is a challenging task because it requires long ded¬ 
icated observations using the instruments with high effective area 
and spectral resolution above 20-30 keV. Only in a fraction of ac¬ 
creting pulsars CRSFs are reliably measured l |Caballero & Wilms| 
|2012||Revnivtsev & Mereghetti|2014) . 

Recently, |Klochkov et al.| ( [2011) and |Reig & Nespoll|(2013) 
pointed out that in addition to the CRSF behaviour, different types 
of the spectrum-luminosity dependence in accreting pulsars can 
be probed by studying the variability of the X-ray continuum with 
luminosity. Specifically, it was found that the X-ray continuum de¬ 
scribed by a power-law dependence F x ~ E r becomes harder with 
increasing X-ray flux in X-ray pulsars showing a positive correla¬ 
tion of the CRSF energy E c with flux, whereas it becomes softer 
or virtually does not change with increasing flux in pulsars show¬ 
ing a negative £’ c (Fx)-correlation (Figs. 3 and 6 in |Klochkov et al.| 
l |2011) ). Clearly, the continuum studies are much less demanding 
from the point of view of data analysis than the spectroscopic anal¬ 
ysis of the CRSF behaviour (if the CRSF is reliably identified in 
the X-ray pulsar spectrum at all). 

The purpose of the present paper is to perform model 2D cal¬ 
culations of the accretion column structure in the radiation diffusion 
approximation and to calculate the emerging X-ray continuum in 
order to see how the hardness of the sidewall emission from accre¬ 
tion columns with different geometries (a filled or hollow cylinder) 
changes with accretion rate variations in the regime of radiation- 
dominated accretion flow braking. We calculate the sidewall en¬ 
ergy flux escaping from the optically thick accretion column and 
the saturated Compton spectrum of extraordinary photons. We find 
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that for typical parameters of X-ray pulsars (the neutron star sur¬ 
face magnetic field B = 3 x 10 12 G, mass accretion rate onto one 
pole M = (1 - 12)x 10 17 g s _1 ), with increasing mass accretion rate 
the hardness of the emergent sidewall emission from the column al¬ 
ways increases. The sidewall emission from the column is Doppler- 
boosted by unbraked matter in the external part of the column to¬ 
ward the neutron star surface. The reflected radiation from the neu¬ 
tron star atmosphere calculated in the single-scattering approxima¬ 
tion is harder than the incident emission, and its addition to the total 
flux form the column renders the spectrum even harder. However, 
starting from some accretion rate the height of the column increases 
such that the fraction of the reflected radiation in the total flux starts 
decreasing. This purely geometrical effect leads to the saturation of 
the hardness ratio at M cr ~ (6 - 8) x 10 17 g s _1 (for the assumed 
neutron star parameters). In the case of the hollow-cylinder accre¬ 
tion column geometry, the saturation of the spectral hardness has 
not been reached in the calculated range of mass accretion rates 
due to smaller photon diffusion times across the columns thickness 
and correspondingly smaller heights of such columns, when the 
contribution of the emission reflected from the neutron star surface 
always dominates. 

The structure of the paper is as follows. In Section[2]we briefly 
describe the observed correlations of CRSF with X-ray flux and X- 
ray continuum properties with X-ray flux. In Section [3] the basic 
equations of the structure of the radiation-supported columns and 
their numerical solution are described. The results of calculations 
of the spectrum are presented in Section [4] which is followed by 
the summary and conclusions. 


2 OBSERVATIONS 


There is a growing number of X-ray pulsars showing the depen¬ 
dence of E c of their X-ray flux (e.g., |Becker et al.] l |2012^ and ref¬ 
erences therein). The cyclotron energy E c in the bright transient 
X-ray pulsar V 0332+53 decreases with luminosity L x (|Tsygankov| 
|et al.|2006||Tsygankov, Lutovinov & Serber|201 Of . A similar nega¬ 
tive £ t .(Lx)-correlation has been reported for another bright tran¬ 
sient pulsar 4U 0115+63 jTsygankov et al.|[2007[ |Mihara, Mak-| 
|ishima & Nagase 2004}. Recently, however, this correlation has 
been questioned by Muller et al.| < [2013| >, who argued that the 
£ c (Lx)- c °rrelation in 4U 0115+63 is an artefact related to the con¬ 
tinuum modelling (see also |Boldin, Tsygankov & Lutovinov| i )20 1 3f 
for an independent analysis). 

There is another type of X-ray pulsars with lower X-ray lu¬ 
minosities and somewhat higher magnetic fields which exhibit a 
positive correlation of E c with Lx ' Her X-l (|Stau bert et a l.|2007| >, 
GX304-1 ( |Yamamo to et al.|201 l| |Klochkov et al,|20 1 2|, andprob- 
ably Vela X-l ((Furst et al.|2014[ and A 0535+26 (| Klochkov et al.| 
[20TT]|Muller et al.|2013) 


As mentioned in the Introduction, a similar bimodality is ob¬ 
served in the dependence of the X-ray continuum hardness on lu¬ 
minosity. This dependence can be studied using, for example, data 
from all-sky monitors such as RXTE/ASM and MAXI, i.e. with¬ 
out dedicated spectroscopic observations which are necessary to 
measure CRSFs. Following this approach, we measured the hard¬ 
ness ratio as a function of luminosity in the accreting pulsars 
GX 304-1, 4U 0115+63, V 0332+53, EXO 2030+375, A0535+26 
and MXB 0656072 using the data from different energy bands of 
RXTE/ASM. The result is shown in Fig.|T] To convert the ASM 
count rates into X-ray luminosities, we used published distances to 
the sources (GX 304-1: ~2 kpc, [Parkes, Murdin & Mason||1980| 


4U0115+63: ~7kpc, |Negueruela & Okazaki||2001| V0332+53: 
~7 kpc, |Negueruela et al|l999| EXO2030+375: ~7kpc, |Wilson| 
|et al.|2002| A 0535+26: ~ 2 kpc, |Steele et al,|1998| MXB 0656- 
072: ~ 3.9 kpc, |McBride et al.||2006 i and their broadband X- 
ray spectra from the archival pointed observations with INTE¬ 
GRAL and RXTE available to us. One can see that at lower fluxes 
the hardness ratio increases with flux. At a certain flux, how¬ 
ever, a flattening of the hardness ratio is observed in 4U 0115+63, 
V 0332+53, EXO 2030+375, A 0535+26 and MXB 0656072. This 
‘turnover’ occurs at the flux roughly corresponding luminosities 
(3 - 7) x 10 37 erg s~'. In GX 304-1, which does not show such 
a turnover, this ‘critical’ luminosity has simply not been reached 
during the outbursts registered by RXTE ASM. 

A similar dependence of spectral properties on the X-ray 
flux at luminosities below ~ 10 37 erg s _1 was also reported from 
spectroscopic observations of the transient Be-X-ray pulsar GRO 
J1008-57 |Ktihnel et al. 1(2014) . 


At the lowest energies covered by RXTE/ASM, the spectrum 
of accreting pulsars is significantly affected by photo-electric ab¬ 
sorption caused mostly by matter in the vicinity of the source. Such 
an absorption naturally influences the derived hardness ratios and 
cannot be accounted for without a high-resolution spectral analysis 
using pointed observations. Neutral absorbing matter surrounding 
the accretor might correlate with the mass transfer rate in the binary 
and might therefore be proportional to the X-ray luminosity of the 
pulsar. Such an effect might at least contribute to the observed pos¬ 
itive correlation between the hardness ratio and flux. However, we 
do not expect any break in the slope of the correlation between 
the absorption column density and flux at L x ~ 10 37 erg s~ ! (al¬ 
though we cannot exclude this possibility). We therefore identify 
the observed breaks with those reported by [Reig & Nespolij ( |20 1 3| > 
at similar X-ray luminosities based on the data from RXTE/PCA 
observations taken at higher energies which should be unaffected 
by the photo-electric absorption. 


3 ACCRETION COLUMN STRUCTURE IN THE 
RADIATION DIFFUSION LIMIT 

3.1 Qualitative considerations 

At high luminosities, L x > 10 37 erg s -1 , the radiation pressure 
dominates, and the braking of the accreting matter flow is due to 
interaction with photons. In this case, the kinetic energy of falling 
particles mpCjl, where v 0 = 10 10 cm s~ ] is the free-fall veloc¬ 
ity near the NS surface, decreases by almost 99% at optical depths 
r ~ (2 — 3)(c/vo) ~ 4 — 9. This was shown for the first time by 
|Davidson 11973J), and in a wide range of parameters was calculated 
by Wang & Frank|(l981 | by solving two-dimensional gas-dynamic 
equations with diffusive radiation transfer. As a result, the shape of 
the zone of energy release looks like a ’mound’ with the character¬ 
istic height zo ~ r 0 . 

With further increasing X-ray luminosity above L x > L*, as 
was first shown by |Basko & Sunyaev] l |1976fr , the height of the 
braking zone starts increasing, and instead of a ’mound’ an opti¬ 
cally thick accretion column appears with zo > r 0 . Most of the ki¬ 
netic energy of the accreting flow is transformed into heat at around 
Zo (which sometimes is referred to as the height of the radiation- 
dominated shock), and below this height the matter settles down 
towards the column base. The thermal energy is adverted by the 
settling matter with an effective velocity of v ~ 0.1 v 0 , and the gen¬ 
erated heat finally escapes through the column sidewalls. 
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It is possible to estimate the critical luminosity L* by equat¬ 
ing the dynamical time of the fall of matter to the diffusion time 
of photons across the column. The characteristic settling time is 
clearly t s = z 0 /v. The characteristic time of the radiation diffusion 
across the filled cylinder column is 

r 7 r 2 0 x ± rgX ± S 

td- — - -P=- (1) 

V d C C V 

where we have eliminated the density p using the mass continuity 
equation 


5 = pv = —r . (2) 

Itr o 

In the case of a hollow cylinder with wall thickness br 0 the diffusion 
time decreases with b: 

b 2 rl b 2 t^x ± b 2 r 2 x ± S bx ± M 

t d - - =s -p =-=- (3) 

Vd C C V C JIV 

By equating t s = t d and noting that the velocity is crossed out, we 
arrive at the expression 

Zo.f =—M, Zo,h = —M (4) 

JtC JtC 

for the filled and hollow cylinder geometries, respectively. 

Now, by equating zo = to, we find the critical mass accretion 
rate M * and the corresponding X-ray luminosity 


V = O.lAfY = 


bx 


-O.lcVo « 2.36xl0 36 [erg s -1 ] ^—] 

V I Oxm/ 


( r ° ) 

1 bx ± ] 

110 5 cm/ 

Ur] 


Here the factor b should be set one for the filled cylinder geometry. 

It is important to realize that the accurate value of L* should 
be found from two-dimensional numerical calculations of the radi¬ 
ation transfer equations, and, for example, |Wang & Frank| ( | 1981 ^ 
found 


V as 5.25 x 10 36 [ergs~‘] ^ 



( 6 ) 


Here the factor (EJE) takes into account the effective change in 
the photon scattering cross-section in the strong magnetic field and 
the deviation of the magnetic field line geometry from the simple 
cylindrical case near the NS surface. This factor is — 1 if the mean 
photon energy E > E c . It is seen that our simple estimate Eq. |5| 
derived above is consistent with this value to within a factor of 
two. Recently, the critical luminosity was reassessed by |Mushtukov| 
|et al. | i ( 20T5) who found it to vary by almost an order of magnitude 
(depending on the magnetic field) within an interval centred at ~ 
10 37 erg s -1 . 

Therefore, the qualitative considerations presented above sug¬ 
gest the increase in the accretion structure height zo x M above 
some critical luminosity ~ 10 37 erg s _1 (see also a more detailed 
derivation in [Arons] <] 1 992ft. Section 3.). 


3.2 Numerical simulations of axially symmetric accretion 
column 

Qualitative considerations described above can be made more pre¬ 
cise by numerical calculations of steady-state accretion columns. 
This radiation-hydrodynamic problem, complicated by the need to 
calculate the radiation transfer in a strong magnetic field, has not 
been solved self-consistently as yet. However, the structure of an 
optically thick accretion column at luminosities above 10 37 erg s _1 
can be calculated in the radiation diffusion approximation with 


grey scattering coefficients along and perpendicular to the mag¬ 
netic field lines i )Davidson]| 1973 [ |Basko & Sunyaevj|1976| |Wang| 
|& Frank|1981| >. The emergent radiation spectrum will be formed 
due to scattering of photons on electrons in the optically thin outer 
layers, which is a separate problem (see the next Section). 

We shall consider two possible axially symmetric geometries 
of accretion columns: a filled cylinder of radius r 0 and a hollow 
cylinder of the inner radius r 0 and thickness br 0 . In both cases r 0 is 
determined by the magnetosphere radius Ra, and in the second case 
the wall thickness br 0 can be related to the distance the disc matter 
enters the magnetosphere before being completely frozen into the 
magnetic field. We shall assume a purely dipole magnetic field and 

the disc matter frozen at a fixed fraction of the Alfven radius. 

. 2 

From the definition of the Alfven radius, R A oc 7, and the 
expression for the polar cap radius ro we see that the polar cap 
radius scales with accretion rate as r 0 oc Mi. The freezing depth 
of disc plasma at the Alfven radius is fixed as A = O.IR 4 , and 
since ^ = b we have b - 0.1. This means that with changing 
mass accretion rate the thickness of the accretion wall in the case 
of hollow cylinder geometry changes as br 0 oc M 1 [j] 

Note that the value of b is difficult to calculate precisely; it 
is clear that it should depend on many factors, including the NS 
magnetic field axis misalignment, the details of the matter entering 
the magnetosphere due to various instabilities, etc., which are be¬ 
yond the scope of the present study. However, we should note that 
in the case of thin cylinder wall, b <K 1 , the radiation diffusion time 
Eq. 0 is proportional t d oc b, thus decreasing the height of the col- 
.(5) 

umn zo at a given mass accretion rate M. This effect is confirmed 
by our calculations (see below). The neutron star magnetic field is 
assumed to be homogeneous across the column. This approxima¬ 
tion is justified for not very high accretion heights (especially in the 
case of the hollow cylinder geometry). 


3.2.1 Boundary conditions 

We are working in cylindrical coordinates r,tp,z centred at the 
columns axis and z = 0 at the neutron star surface. The initial veloc¬ 
ity of the falling matter at high altitude above the surface of the neu¬ 
tron star is Vo = 10 10 cm s _1 . At the cylinder base (z = 0), the ve¬ 
locity is v = 0. As a boundary condition at the column side surface 
we used the relation between the radial energy flux to the energy 
density U in the form F,.(r 0 ,z ) = 2cU(r 0 ,z)/3, which roughly cor¬ 
responds to the conditions expected in the scattering atmospheres 
in the Eddington approximation. 


3.2.2 Basic equations 

The cylindrical symmetry of the problem makes it essentially two- 
dimensional. The steady-state momentum equation (ignoring grav¬ 
ity, which is very good approximation as discussed e. g. in|Basko &| 
|Sunyaev[ l jl976| l) for accretion braking by the radiation with energy 
density U reads 

(S-V)v = ~VI/, (7) 

where S = p\ = const is the mass continuity equation. The integra¬ 
tion of these equation yields 


1 In principle, one may consider further complications, e.g., using a specific 
model of the magnetosphere-disc coupling, like it was done in|Mushtukov| 


et al. 


2015 


This only slightly changes the dependence of bro on M. 
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a) b) 




Figure 2. Contours of constant Q (0.9 to 0.1, from top to bottom). Top panels: the hollow cylinder geometry with b = 0.1 calculated for Myj — 1 (a) and 
M\n = 5 (b). Bottom panels: the filled cylinder geometry for M\-j = 1 (c) and M\-j = 5 (d). 


(7 = 35(v 0 - v). (8) 

Following |Davidson| ( |1973) l, the energy equation can be writ¬ 
ten as 

V • F = -S • V j, (9) 

where we has neglected the flow of internal energy of the falling 
matter. 

The radiative transfer equation in the diffusion approximation 
is 

F = ——VLC + -U\, (10) 

3 xp 3 

where the inflow of gas is considered to be stationary. 


From equation in cylindrical coordinates we find the com¬ 
ponents of the radiation energy flux: 


c dU 
3x±p dr 


1 + 


1 1 \dU\ 

3x±p U I dr I 


f 7 = 


_ _ 4 jj 

3x,|p dz 3 U 

i . 1 1 \dU\ ' 

1 3x l]P U I 8z I 


(ID 


( 12 ) 


where the coefficients (l + |f |) ‘ and (l + |f |) 

are introduced to use the modified diffusion approximation in the 
optically thin regions where the change of the radiation energy den¬ 
sity on the photon mean free path becomes significant (see |Wang] 
|& Frank] t |1981| >). Equations {7} - i | 1 2| i describe the accretion 
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shock mediated by radiation in the diffusion limit and provide the 
structure of the zone where most of the accretion energy respon¬ 
sible for the formation of the emergent spectrum is released. The 
’sinking zone’ downstream the radiation shock, where the matter 
settles down subsonically and advects thermal energy towards the 
NS surface, should be treated separately (see e.g. |Basko & Sun-| 
|yaev| p976 I). Different physical effects can appear in this zone, e.g. 
’photon bubble oscillations’ discussed in Arons (19921; Klein et al. 
(1996). 


3.2.3 Equations in dimensionless form 


The convenient variables are dimensionless velocity of matter 
squared Q = v 2 /Vq, the dimensionless column radius r = 
x T Sor/c = T T vr/c and height z = x t Sqz/c = r T vz/c , where 
S o = S(Mn = 1) = 8.79 x 10 6 ergs -1 cm -2 (here and below 
Mn = M/(10 17 gs -1 ) is dimensionless mass accretion rate onto 
one pole). For the assumed value of rg(Mn = 1) the dimension¬ 
less factors almost exactly equal to 10 -4 cm -1 . For the scattering 
cross-sections we shall assume x ± = x T and Xj/x y = 10. 

With this notation and dimensionless variables, the system of 
equations & & {HI- (12) is equivalent to one non-linear partial 
differential equation (cf.|Davidson|(l97~3)) 


K ± 


d 2 Q 1 dQ 
dr 2 + f dr 


+ K,| 


x T d 2 Q 
x\\ dz 2 


8 ag(VG-e) 


dQ 

- # = 0, (13) 
dz 


where the coefficients K ± = (l + |ff|) and ^11 = 

11 + — - ,, vo — 1^1) are calculated at each consecutive iteration 

\ *n 6c(i-V5) i & 7 

using values of Q obtained in the previous iteration. 

Equation (13) was solved numerically using the finite differ¬ 
ence approximation method. A time-relaxation scheme was used. 
In this approach, the equation in the form L[Q} = 0, where L 
is the corresponding differential operator, has been replaced by a 
parabolic equation of the form ^ - L[Q] = 0. The steady-state so¬ 
lution of the last equation coincides with the solution of the original 
elliptic equation. An explicit finite-difference scheme was used to 
solve the parabolic equation. 


3.3 Results of numerical calculations 

Numerical simulations were carried out for accretion rates onto one 
pole in the range 10 17 - 1.2 x 10 18 g s -1 for the hollow cylinder 
geometry and the filled cylinder geometry. The column structure at 
different accretion rates is shown in Fig. [2] Note that at a given mass 
accretion rates and the NS magnetic field, the height of the braking 
zone in the hollow cylinder case is smaller than that in the case of 
the filled cylinder, as expected from the dependence of the radiation 
diffusion time on the wall thickness oc b. Clearly, in very geometri¬ 
cally thin walls, b/r 0 <K 1, the flow will be optically thin as well (at 
least for the extraordinary mode with angle-independent scattering 
cross-section in the strong magnetic field). This may suggest that 
even if the flow was initially confined within a thin layer, near the 
NS surface the huge energy release may change the character of the 
flow. This possibility is worth investigating further. 

Guided by the velocity braking in a radiation-dominated 
shock, v « l/7v 0 , we define the ’height’ of the accretion mound 
Zo(M) as the distance at the centre of the structure where the gas 
flow velocity is reduced to one seventh of the initial value. The 
(dimensionless) height calculated for different accretion rates were 
fitted by a power law 



Figure 3. The radial energy flux F, (t = 2/3, z) calculated along the height 
of a filled cylinder column at different mass accretion rates. Half of this flux 
(extraordinary photons only) escapes sidewall. The lines from bottom to up 
correspond to M\i = 2. 3, 5, 8 and 12, respectively. 


04) 

where z 0 = x T SoZo/c. We found that in both geometries a ~ 1 and 
/3 > 0 (a - 1.1, (3 - 0.3 for the case of hollow cylinder and a - 0.9, 
/( - 13.5 for the filled cylinder). Therefore, the numerical solution 
behaves in agreement with the qualitative physical considerations 
about the structure of optically thick accretion mounds described 
above. 


4 CHANGE IN THE HARDNESS RATIO FROM 
ACCRETION COLUMNS 


4.1 Saturated Compton spectrum of sidewall emission 

The LTE treatment of the spectrum of the sidewall emission from 
accretion columns is clearly a strong oversimplification. To take 
into account the scattering on electrons in the optically thin bound¬ 
ary of the column, we can use the model calculations of the radi¬ 
ation transfer problem in a semi-infinite plane-parallel atmosphere 
with strong magnetic field. In this case the emerging spectrum will 
be formed in the saturated Compton regime, with the mode 1 pho¬ 
tons (extraordinary, i.e. polarized perpendicular to the k - B-plane, 
where k is the photon wavevector and B is the magnetic field vec¬ 
tor) predominantly escaping through sidewalls (Lyubarskii| 1986) . 
The intensity of mode 2 (ordinary, i.e. polarized in the k-B plane) 
photons is comparable to that of mode 1 photons only at large an¬ 
gles to the normal, // ~ 0, and thus insignificantly contribute to the 
total flux from the unit surface area due to the geometrical factor. 

In this regime, the specific intensity of extraordinary photons 
at energies far below the cyclotron resonance can be written to an 
accuracy of a few per cents as (Lyubarskii|1986) 

2V3 , hvv 2 hv 

I = (15) 

5 c-T 0 

where // = cos 6' is the angle between the normal to the atmosphere 
and the escaping radiation direction in the plasma reference frame. 
To is the characteristic optical depth of the problem, v g = 2S 
the electron gyrofrequency in the magnetic field. By expressing the 
characteristic optical depth r 0 through the emergent radiation flux 
using Eqs. (34-37) from |Lyubarskii[fT986) , we find 




hvY O 


— e kT . 

10jt \kT I v 


(16) 
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Figure 4. Left: The spectrum of sidewall emission from optically thick filled cylinder accretion column Jl 6 j for mass accretion rates M \7 = 2, 3, 5, 8 and 12 
(from bottom to up, respectively). Right: The spectral hardness ratio HR as a function of mass accretion rate. 


Therefore, we are in the position to calculate the emerging 
spectrum from the column in this approximation using the solutions 
for the column structure obtained above. To do this, it is sufficient to 
substitute the total radial energy flux at each height of the column 
<!> —> F r (z)/ 2, and temperature T —> T(z ) estimated deep inside 
the column (at the optical depth ~ 1), so that T(z ) = (I/(z)/a,.) 1/4 ). 
Note here that at large optical depth the number density of ordi¬ 
nary and extraordinary photons are equal, with the extraordinary 
photons becoming overwhelmingly dominant in the optically thin 
outer layers. Therefore, only half of the radial energy flux, F r (z)/ 2, 
should be taken into account in the spectral calculations. The sec¬ 
ond half of the radial flux (in ordinary photons) is admixed to the 
vertical energy flux component, F z . Some fraction of this energy 
flux can propagate (mostly as ordinary photons) along the mag¬ 
netic field into the optically thin upper part of the column; we do 
not calculate their fate in this paper and restrict ourselves by this 
quantitative note. The radial energy flux F r (z) from the filled col¬ 
umn is shown in Fig.[3]for different mass accretion rates. 

The resulting spectra of the sidewall emission from the opti¬ 
cally thick accretion column integrated over the column height z is 
shown in Fig. [4] (left panel) for the case of the filled column and 
different M. The spectral hardness ratio 

12 keV 

f Cctv 

5 keV 

HR=-— -, (17) 

3 keV 

f F v dv 

1.3 keV 

is shown in the right panel of this figure as a function of M. It is seen 
that the hardness ratio of the spectrum of the column monotonically 
increases with mass accretion rate. 


4.2 Account for the reflected component 

We have seen that the saturated-Comptonized spectrum of sidewall 
emission from optically thick accretion column gets harder with 
increasing M, contrary to what is observed and discussed above. 
However, an important feature of accretion in X-ray pulsars should 
be added at this point. Namely, we should take into account the 
fact that the electrons in the optically thin part of the column are 
moving with high velocity Vo ~ 1/3c, and thus the emission beam 



Figure 5. The single-scattering Compton X-ray albedo from a plane- 
parallel neutron star atmosphere with strong magnetic field. The electron 
number density is n e = 5 X 10 25 cm' 3 , electron temperature T = 3 keV, 
magnetic field B = 3 X 10 12 G, photon incident angle to the magnetic field 
ji/ 4. Extraordinary and ordinary photon spectra are shown by the solid and 
dashed lines, respectively. The cyclotron line is ignored. 

of the column 0 should be Doppler-boosted towards the neu¬ 
tron star surface and reflected from the neutron star atmosphere, 
as discussed in |Lyubarskii & Syunyaev[ ( jl988| l and Poutanen et al. 
l |2013| l. The reflection coefficient in the case of single Compton 
electron scattering in strong magnetic field is determined by the ra¬ 
tio A(v) = k sc /(k sc + K abs ), where k sc and K ab s is the scattering and 
absorption coefficients for photons in the strong magnetic field, re¬ 
spectively. The reflection coefficient (X-ray albedo) is calculated in 
the Appendix |A| and is presented in Fig. [5] Note that the calculated 
X-ray albedo turns out to be very similar for both extraordinary and 
ordinary photons and is virtually insensitive to the photon incident 
angle to the magnetic field. Everywhere below we shall use X-ray 
albedo for extraordinary photons only. It is seen that soft photons 
are mostly absorbed, while hard photons are reflected by the neu¬ 
tron star atmosphere. Therefore, the reflected spectrum is harder 
than the incident one. 

Let a be the angle between the photon wavevector and the 
plasma bulk velocity vector. Due to the boosting, part of radiation 
from the column will be intercepted by the neutron star surface and 
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reflected. The critical angle or* within which the radiation from the 
height z above the surface will be intercepted by the neutron star 
(in the Schwarzschild metric) is (|Poutanen et al. 2013} 


RnS j 1 ~ r s/(RNS + z) 
RnS + £ \ l _ rs/^AtS 


(18) 


where r$ = 2 GM/cr is the Schwarzschild radius. The difference be¬ 
tween the polar axis and the column cylinder side can be neglected 
for r 0 « R ns . 

Let £ be the surface of the column, £'(</>, z) be a point on the 
surface and rf£ = r 0 difidz be the elementary area. The flux from the 
surface element at frequency v (in conventional notations) reads 


/,(£') = 


dE(L') 

dtdvrodtfidz 


(19) 


and can be separated into three parts in accordance to the angle 
relative to the observer. Photons escaping with angles a > a r* to the 
column are directly seen as the proper column radiation f" jl . The 
second part includes the radiation f" 1 intercepted by the neutron 
star and reflected from the neutron star surface with X-ray albedo 
T(v). The third part f° bs is the radiation absorbed and re-radiated 
by the neutron star atmosphere. Then we can write: 


flux. This, in turn, depends on the height of the column. We have 
calculated the total flux as the sum of the direct and reflected com¬ 
ponent as a function of M with taking into account the photon ray 
propagation in the Schwarzschild metric of the neutron star with a 
fiducial mass M = 1.5M 0 and radius R = 10 and 13 km. The mag¬ 
netic field of the neutron star is set to 3 x 10 12 G, so that the CRSF 
energy is about 35 keV. The result is presented in Fig. [6] It is seen 
that the hardness ratio HR first increases with the mass accretion 
rate (X-ray luminosity), but starting from M cr ~ (6 - 8) x 10 17 g s -1 
it gets saturated and even decreases. This is in agreement with ob¬ 
servations shown in Fig. |T] 

For the hollow cylinder accretion columns which have smaller 
height (see Fig. [2}, the fraction of the reflected component in the 
total emission does not virtually change with mass accretion rate 
within the calculated range, therefore the hardness ratio of the total 
spectrum monotonically increases with M in this range (see Fig.[7]l. 
Note also that in this case the spectrum is harder than in the case of 
the filled column because the reflected (harder) emission dominates 
in the total spectrum. 


5 SUMMARY AND CONCLUSIONS 


/,,(£') = f d P 1 ‘ - d cos a 

J 1 7 J a cos a 

-1 

cos a* 1 

_ r dfv(L r dfvd. ) d cos a 

J a cos a J a cos a 

— 1 cos a* 

cos a* 

C dfy(L') 1 

= -f^-d cos a 

J a cos a 
-1 

1 1 
+/l(v) f d / v(L ’ d cos a + (1 — /l(v)) f d / y(1, ^ d cos a 

v 7 J d cos ce v v 7/ J d cos a 

cos O'* cos a* 

= /;■"(£') + /; e/ (£') + /; is (£'). 


( 20 ) 


The integrands (Poutanen e t al.|2013} in the cylindrical coordinates 
with account for the axial symmetry (7 V (£') = 7 v (z), giving by C3) 
have the form: 


dfW) 

dcosa 


2D 3 I nD 
= Iv(z) -sino-^l H—— sin n 


( 21 ) 


where the Doppler factor is D = 1 /(y(l - p cos or)), y = 1/ y/l - ft 2 
is the plasma Lorentz factor and /? = v/c. Here we have taken into 
account that the specific flux (per unit frequency interval) is D times 
as small as the integral flux. Then for the total radiation from the 
column at frequency v we obtain: 


U = ff fy(P)dL = 2jtr 0 JfAz)dz 

£ 0 

= 2itr 0 J\r'(z) + f7 f (z) + f v abs (z))dz 

o 

= Lf + L[ ef + L abs , (22) 

where z m ax is the upper limit of the computational area. The ob¬ 
served X-ray flux from one column is F v = ( L c v ° 1 + L r v e ^)/(4nd 2 ), 
where d is the distance to the source. For the illustrative purposes, 
we assume <7 = 5 kpc. Taking into account of the radiation absorbed 
by the neutron star atmosphere L abs will add more soft photons to 
the total spectrum. 

Since we are observing both the direct and reflected radiation 
from the column, the total change in the hardness of the spectrum 
will depend on the fraction of the reflected radiation in the total 


In this paper, using RXTE/ASM archival data we studied the be¬ 
haviour of the spectral hardness ratio as a function of X-ray lumi¬ 
nosity at different accretion rates in a sample of six transient X-ray 
pulsars (EXO 2030+375, GX 304-1, 4U 0115+63, V 0332+63, A 
0535+26 and MXB 0656-072). In all cases, the hardness ratio is 
found to increase with X-ray flux at low luminosities and then to 
saturate or even slightly to decrease above some critical luminosity 
of a few times 10 37 erg s~' This behaviour confirms earlier findings 
by |Reig & Nespoli|j20 1 3[ >, and |Klochkov et al7] < [201 1| > and corre¬ 
lates with the behaviour of the CRSF energy with flux in the CRSF 
sources. 

At low mass accretion rates and low X-ray luminosities, the 
braking of infalling matter is expected to be mediated by Coulomb 
interactions (|Zel’ dovich & Shaku raj |1969| |Nelson, Salpeter~&| 
|Wasserman| 1993; Bec ker et al.|2012) and occurs close to the neu¬ 
tron star surface, essentially at a height of the homogeneous at¬ 
mosphere above the neutron star surface, without the formation of 
an accretion column. When the X-ray luminosity increases above 
some critical value ~ 10 37 erg s 1 i DavidsonJl973 Basko & Sun- 


|yaev|1976[|Wang & Frank|1981||Mushtukov et al. 2015) , the radia- 
tion pressure of generated photons decelerates the infalling matter. 
An optically thick accretion column arises with the height increas¬ 
ing almost linearly with accretion rate. 

Therefore, at low luminosities < L 
emerging X-ray spectrum is essentially formed by ordinary pho¬ 
tons in a magnetized optically thin atmosphere. The hardness of 
the spectrum is determined by the Comptonization y-parameter, 
y = (kT/m e c 2 )T, which increases with accretion rate. This implies 
that in this case the hardness of the emerging radiation will increase 
with luminosity. 

At higher luminosities an optically thick accretion column 
above the magnetic pole arises, and thermal radiation generated 
deep inside the column escapes sidewall predominantly as extraor¬ 
dinary photons. Since at energies far less than the CRSF energy 
the electron scattering cross-section strongly decreases with fre¬ 
quency, the spectrum of the emergent continuum radiation will be 
softer than the diluted Wien spectrum expected to be formed in the 
saturated Compton regime. 


10 37 erg s 1 the 


-hv/kT 


Lyubarskii (19861. 


We have calculated the structure of axially symmetric opti- 
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Figure 6. Left: the total spectrum of direct sidewall and reflected from the neutron star atmosphere from the optically thick filled accretion column for mass 
accretion rates Myj = 2, 3, 5, 8 and 12 (from bottom to up, respectively). Right: the hardness ratio HR of the total spectrum. Shown are calculations for the 
neutron star radius R^s = 10 km (squares) and R^s = 13 km (circles). 
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Figure 7. Hardness ratio of the total emission (direct plus reflected from 
the neutron star atmosphere) from a hollow cylinder accretion column with 
b = 0.1 /■(} as a function of mass accretion rate M. 


cally thick accretion columns in the diffusion approximation for 
radiation transfer for different mass accretion rates. Two accretion 
column geometries, filled cylinder and hollow cylinder, were con¬ 
sidered. In both cases we have found the expected almost linear 
growth of the column height with mass accretion rate. The con¬ 
tinuum X-ray spectrum of the sidewall emission in both geometries 
gets harder with increasing mass accretion rate. Taking into account 
of the radiation reflected from the neutron star atmosphere allowed 
us to reproduce the observed saturation of the continuum hardness. 
This is purely geometrical effect: with increasing mass accretion 
rate the height of the column increases, so does the fraction of the 
reflected radiation. However, starting from some height, the frac¬ 
tion of the reflected radiation in the total flux stops increasing and 
starts decreasing. Since the hardness of the reflected radiation is 
higher than that of the incident emission (because of the strong en¬ 
ergy dependence of the scattering and absorption coefficients in the 
strong magnetic fields), the hardness of the total spectrum (direct 
sidewall emission from the accretion column plus reflected radi¬ 
ation from the neutron star atmosphere) saturates (and even can 
slightly decrease) beyond some characteristic mass accretion rate 


~ (6 - 8) x 10 17 erg s _1 for the typical neutron star parameters: 
M ns = 1.5M 0 , R ns = 10 km and 13 km, B = 3 x 10 12 G. 

In the case of the hollow cylinder accretion column geometry 
with a fiducial wall thickness of 0.1 the outer column radius r 0 , we 
have not found the spectral hardness ratio saturation with increas¬ 
ing mass accretion rate up to the maximum value of our calcula¬ 
tions 1.2 x 10 18 g s -1 . Apparently, the height of the column in this 
case is too low, the hard reflected component dominates at all mass 
accretion rates, and the (softer) direct sidewall emission remains 
subdominant in the total spectrum. 

Therefore, our model calculations lead to the following con¬ 
clusions. 

1. The spectral hardening in X-ray pulsars with positive CRSF 
energy dependence on X-ray flux can be explained by increasing 
of the Comptonization parameter y in the slab atmosphere of the 
accretion mound. It is in this regime that positive correlation of the 
cyclotron line is observed, for example, in Her X-l ([Staubert et al.| 
[2U07l l. 

2. At high accretion rates, the radiation-supported optically 
thick accretion column grows above the polar cap, the sidewall 
emission from the column is formed by extraordinary photons in 
the saturated Compton regime. The spectrum of this emission gets 
harder with increasing mass accretion rate. 

3. With further increasing mass accretion rate the height of the 
column increases, such that the fraction of radiation reflected from 
the neutron star atmosphere starts decreasing. As the reflected ra¬ 
diation is harder than the incident one, the spectrum of the total 
emission (direct plus reflected) stops hardening (and even becomes 
slightly softer). This happens at the mass accretion rate onto one 
pole M cr ~ (6 - 8) x 10 17 g s _1 , in qualitative agreement with ob¬ 
servations (see Fig.[IJ. 

4. In the frame of this model, the saturation of the spectral 
hardness in the case of a hollow cylinder geometry of the accretion 
column can be achieved at much higher accretion rates (roughly, 
scaled with the relative thickness of the column wall, r 0 /b), because 
the characteristic height of the column in this case is correspond¬ 
ingly smaller than that of the filled column. 

5. Variations in the pulse profile form and polarization 
properties during the transition from the Coulomb-braking low- 
luminosity regime to the radiation-braking high-luminosity state 
should occur with changing mass accretion rate (X-ray luminosity) 
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in accreting X-ray pulsars. The profile change has already been re¬ 
ported during the low-states in Vela X-l ^Doroshenko, Santangelo| 
|& Suleimanov|20TT| >. 

We conclude that X-ray continuum observations in transient 
X-ray pulsars can be used as additional diagnostic of different ac¬ 
cretion regimes and column geometry near the surface of highly 
magnetized neutron stars. 
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APPENDIX A: X-RAY ALBEDO FROM NEUTRON STAR 
ATMOSPHERE 


The X-ray albedo (the fraction of the reflected to incident flux) from 
the neutron star atmosphere in strong magnetic field in the gen¬ 
eral case depends on the angle between the incident photon and the 
magnetic field and the photon polarization. Not solving exactly the 
problem, we calculate single-scattering albedo for a plane-parallel 
atmosphere T J (v) (J = 1,2 for extraordinary and ordinary photon 
polarizations, respectively): 


T'(v) = 


■ + P,. 


(Al) 


where k[, and k 1 . are scattering and absorption coefficients in the 
magnetic field, respectively. 

We calculate the coefficients following|Venturaj (J 1979^; Nagel 
(1980u Ho & Lai| ( [2001| >. The electron and ion scattering coeffi¬ 
cients for mode j reads: 


«ic,e = "eO-T J] 


pi I 2 


^ (at + acj c e ) 2 


(A2) 


Z~m e 

Am,, 


HjCy T 


Z a>~ 

, (a> + aut c ,■)- 

tt=-l ’ 


J I 2 


K\ 


(A3) 


where <r T is the Thomson scattering cross-section, m e and m p is 
the electron and proton masses, n e and n, is electron and ion con¬ 
centrations, at ce and at c j is the electron and ion cyclotron angular 
frequency, e 2 _ v e\, are e J Q are the components of the photon polar¬ 
ization vector e ; for mode j (see subsection 2.5 and Appendix A in 
l|Ho & Lai 1 2001 [i). The absorption coefficients are 


= z 


pi l 2 c 


(at + aat c<e ) 2 


(A4) 


j _ (Z 2 m e \ 2 % n a > 2 j 2 
^ ( Am,, j 7? Zj (w + au c ,) 2 ^ 8 ‘ 

' ^ ' a ——1 

where 


(A5) 
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k 0 = 4n 2 Z 2 a 3 



_ tkj \ 

- e w I 


(A6) 


is the free-free absorption coefficient in the absence of a magnetic 
field, g ± i and go are the modified Gaunt-factors (|Nagel| 1980|), 


g±i = J' exp ^ sinh 2 xj C, | e 2jI J dx, (A7) 

go = J' exp ^ sinh 2 J.J 2^- e 2x C 0 | e 2x j dx. (A8) 


Here a> c is electron (ion) cyclotron angular frequency. Functions Cj 
and Co are calculated in |Virtamo & Jauho] ( |1975| ). The total coeffi¬ 
cients are: 


(A9) 


k j . = k 1 , + k 1 . .. 

abs abs, e abs, i 


(A10) 


The coefficients k sc and K abs depend on many factors. They 
are determined by the characteristic of reflecting atmosphere, in¬ 
cluding the chemical composition, ion number density and temper¬ 
ature, as well as by the photon incident angle to the magnetic field 
and the magnetic field strength. The last two parameters turn out to 
be weakly changing the X-ray albedo A J (v) in the interesting range 
of parameters (see Fig. |5j. In the calculations we have assumed 
purely hydrogen atmosphere with ion charge Z = 1, the mass num¬ 
ber A = 1 and the ion number density n, = 5 x 10 25 cnT 3 (it can be 
higher in deeper layers ( |Ho & Lai|2001[ l). The electron temperature 
is fixed at T = 3 keV, which corresponds to the expected Compton 
temperature during reflection of X-ray emission form the column. 

The single-scattering X-ray albedo as a function of frequency 
for ordinary (j = 2) and extraordinary (j = 1) photon polarizations 
is shown in Fig.|5] 
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